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Abstract 
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The investigation of the solutions for quantum mechanical systems having certain poten- 
tials in the case of position-dependent mass (PDM) [1, 2] has been received great attentions. 
Many authors have studied the solutions of different potentials for spatially-dependent mass, 
such as hypcrgcometric type potentials [3], Coulomb potential [4], PT-symmetric kink-like, 
and inversely linear plus linear potentials [5]. It is well known that the theory based on 
the effective-mass Schrodinger equation is a useful ground for investigation of some physical 
systems, such as semiconductor heterostructures [6], the impurities in crystals [7-9], and 
electric properties of quantum wells, and quantum dots [10]. In the present work, we tend 
to solve the Dirac-Morse problem within the PDM formalism. 

The pseudospin symmetry is an interesting result appearing in Dirac equation of a particle 
moving in an external scalar, and vector potentials in the case of it when the sum of the po- 
tentials is nearly zero. It was observed that the single particle states have a quasidegeneracy 
labeled with the quantum numbers and s, which are called the pseudo-orbital angular mo- 
mentum, and pseudospin angular momentum quantum numbers, respectively [11-16]. The 
concept of pseudospin symmetry has received great attentions in nuclear theory because of 
being a ground to investigate deformation, and superdeformation in nuclei [17, 18], and to 
build an effective shell-model coupling scheme [19, 20]. The symmetry appears in that case, 
when the magnitude of scalar potential is nearly equal to the magnitude of vector potential 
with opposite sign [14, 21-25] and the Dirac equation has the pseudospin symmetry, when 
the sum of the vector, and scalar potentials is a constant, i.e., E(r) = Vv{r) + Vs{r) — const. 
or dTi{r)/dr — [16]. The spin symmetry is another important symmetry occurring in Dirac 
theory in the presence of external scalar, and vector potentials. The spin symmetry appears 
in the Dirac equation, when the difference of scalar, and vector potentials is a constant, i.e., 
A(r) = V;(r) - Vs{r) = const. [14, 16]. 

Recently, the pseudospin and/or spin symmetry have been studied by many authors for 
some potentials, such as Morse potential [26-28] , Woods-Saxon potential [29] , Coulomb [30] , 
and harmonic potentials [31-33], Eckart potential [34-36], Poschl- Teller potential[37, 38], 
Hulthen potential [39], and Kratzer potential [40]. In Ref. [41], the bound-state solutions 
of Dirac equation are studied for generalized Hulthen potential with spin-orbit quantum 
number k in the position-dependent mass background. In this letter, we tend to show that 
the new scheme of the Nikiforov-Uvarov (NU) method could be used to find the energy 
spectra, and the corresponding eigenspinors within the framework of an approximation to 
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the term proportional to for arbitray spin-orbit quantum number k, i.e. k ^ 0, when 
the mass depends on position. The NU method is a powerful tool to solve of a second order 
differential equation by turning it into a hypergeometric type equation [42]. 

Dirac equation for a spin-^ particle with mass m moving in scalar Vs{r), and vector 
potential Vy (r) can be written as (in ^ = c = 1 unit) 

[a . P + /3(m + Vs{r))] ^„,(r) = [E - V,{r)] ^nn{r) • (1) 

where E is the relativistic energy of the particle, P is three-momentum, a and /3 are 4x4 

/O ct\ / /\ 

Dirac matrices, which have the forms of a = and /3 = { , respectively, [43] . 

\a OJ \-I OJ 

Here, u is a three-vector whose components are Pauli matrices and / denotes the 2x2 unit 
matrix. J denotes the total angular momentum , and K = —f3{a.'L + 1) corresponds to the 
spin-orbit operator of the Dirac particle in a spherically symmetric potential, where L is 
the orbital angular momentum operator of the particle. The eigenvalues of the spin-orbit 
operator K are given as k, — ±{j + l/2), where k, — —{j + 1/2) < correspond to the aligned 
spin j — £ + 1/2, and k— {j + 1/2) > correspond to the unahgned spin j — i — 1/2. The 
total angular momentum quantum number of the particle is described as j = ^ -|- s ,where 
£ = i + 1 is the pseudo-orbital angular momentum quantum number, and § = 1/2 is the 
pscudospin angular momentum quantum number. For a given k, = ±1, ±2, . . ., the relation 
between the spin-orbit quantum number k , and " two" orbital angular momentum quantum 
numbers are given by k,{k, + 1) = i{i + 1), and k,{k, — 1) = i{i + 1). 

The Dirac spinor in spherically symmetric potential can be written in terms of upper and 
lower components as 



Xnn {r)YU9, 0) 
' V i<l>n. {r)Yfj9, <f>) 



^nn{r)=- ""^^'^ T^^' , (2) 



where Yj^{9, 0), and Yj^{9, (j)) are the spherical harmonics, and Xnn {'>^)/'>^, and (pnK {r)/r are 
radial part of the upper and lower components. Substituting Eq. (2) into Eq. (1) enable us 
to write the Dirac equation as a set of two couple differential equations in terms of Xuk (t^) 
and (f)ni^ (r). By eliminating Xuk ('") or (p^K (i^) in these coupled equations, we obtain 

f (P k(k + 1) 1 /dm{r) dA{r)\ d k ^ ,^ / n,^ / x / x /on 

{ - + - ^) + = J^A(r)M.(r)x„(r) ,(3) 
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f (P k(k — 1) 1 /dm(r) dT,(r) \ , d /tx-i, , . , ^ , 



where MA(r) = m + - A(r) , Ms(r) = m - K« + S(r), and A(r) = K (r) - Vs (r), 
S(r)-K(r) + K(r). 

In the NU-mcthod, the Schrodinger equation is transformed by using an appropriate 
coordinate transformation 

(72(s)^"(s) + cT(s)f(s)*'(s) + a(s)^(s) = , (5) 

where (7(s), 5'(s) are polynomials, at most second degree, and f(s) is a first degree polyno- 
mial. The polynomial 7r(s), and the parameter k are required in the method 



As) = \ Wis) - f{s)] [a'{s) - f{s)Y - a{s) + ka{s), (6) 

A = A; + 7r'(s), (7) 

where A is a constant. The function under the square root in the polynomial in 7r(s) in 
Eq. (6) must be square of a polynomial in order that 7r(s) be a first degree polynomial. 
Replacing k into Eq. (6), we define 

r(s) =f(s) + 27r(s). (8) 

where the derivative of r(s) should be negative [42]. Eq. (5) has a particular solution with 
degree n, if A in Eq. (7) satisfies 

A = A„=-„.'-Mri^, „ = 0,1,2,... (9) 

To obtain the solution of Eq. (5) it is assumed that the solution is a product of two 
independent parts as ^(s) = (j){s) y{s), where y{s) can be written as 

where the function p{s) is the weight function, and should satisfy the condition 

[a{s) p{s)]' = t{s) p{s) , (11) 
and the other factor is defined as 

1 (i(/)(s) 7r(s) 



(f){s) ds a{s) 



(12) 



In order to clarify the parametric generalization of the NU method, let us take the following 
general form of a Schrodinger-like equation written for any potential, 

When Eq. (13) is compared with Eq. (5), we obtain 

f{s) = «! - a2S ; a{s) = s(l - ags) ; a{s) = -^is^ + - 6 ■ (14) 
Substituting these into Eq. (6) 



(16) 



tt{s) = 0:4 + 0^55 ± Y (tte ~ /cq;3)s^ + (ay + A;)s + ) (15) 
where the parameter set are 

Q;4 = 1(1 - tti) , "5 = 1 ("2 - Stta) , Q!6 = ai + Ci 
ay = 2ci;4q;5 — , = Q;| + ^3 . 

In NU-method, the function under the square root in Eq. (15) must be the square of a 
polynomial [42] , which gives the following roots of the parameter k 

h,2 = -(ay + Saaag) ± , (17) 

where ag = a^a-j + a\as + cte • We obtain the polynomials 7r(s) and t(s) for k = —{a^ + 
2a3a^) — 2y/a^a9, respectively 

7r(s) = a4 + q;5S - [(i/a^ + a3V^)s - , (18) 

t(s) = ai + 2q;4 - {a2 - 2a^)s - 2 [(v^ + a2,^/a^)s - ^/ai] . (19) 

Thus, we impose the following for satisfying the condition that the derivative of the function 
r(s) should be negative in the method 

t'{s) = -(q;2 - 20:5) - 2(v^ + ag-^/o^) 

= -2a3-2(V^ + a3^/^) < 0. (20) 

Prom Eqs. (7), (8), (19), and (20), and equating Eq. (7) with the condition that A should 
satisfy given by Eq. (9), we find the eigenvalue equation 

Q;2n - (2n + 1)q;5 + {2n ^- l){y/a^ ^- a^^/oi^) ^- n{n — l)a2, 

+ ay + 2Q;3a8 + 2-,/a^ = 0. (21) 



We obtain from Eq. (11) the polynomial p{s) as p{s) — s°'^°~^{l — a^s) "s and 
substituting it into Eq. (10) gives 

yn(.) = - - 2^3.) , (22) 

where aio — ai + 2a4 + 2y/a^, an — a2 — 2ol^ + 2(y^ + OL-iyfoL%) and P^°"^)(l — la-^s) are 
the Jacobi polynomials. Prom Eq. (12), one obtaines 

0(5)^ 5-1.(1 _«3s)-°--^, (23) 

then the general solution ^(s) = 4>{s)y{s) becomes 

^(s) = s"-(l - a3s)""""'^Pi"'°"''^~"'°~'\l - 2a3s). (24) 



where 0:12 = 04 + ^/as and «i3 = 05 — (y^ + as^/og )• Let us study the case where the 
parameter = 0. In this type of problems, the eigenfunctions become 

^(s) = e""' L^'"-\aiis) , (25) 

when the limits lima^-^oPn ' (1 ~ ctS'^) = Ln^°~^i<^iis) ^^^03^0(1 ~ 

a^s) = e"^^* are satisfied and the corresponding energy spectrum is 

— 2a5n + {2n + l){^/aQ — a3^/a8) + n{n — l)a2 + ol-j 

+ 2a2,a^ - 2^/a^ + = Q . (26) 

The generahzed Morse potential is given by [44] 

Vm{x) = L>e-2^=^ - 2De-^'' , (27) 

where x — {r / Tq) — 1 , ^ — avQ , D is the dissociation energy, tq is the equilibrium distance, 
and a is the potential width. The term proportional to 1/r^ in Eq. (4) can be expanded 
about a; = [45] 



^7^= (Y^ = «o(l-2x + 3x^ + ...); ao= 



Instead, we now replace Vm{x) by the potential [45] 

Vm{x) = ao(ai + 026-'^^ + 036"^^^) , (29) 
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Expanding the potential Vm{x) around x — 0, and combining equal powers with Eq. (28), 
one can find the arbitrary constants in the new form of the potential as 

3 3 4 6 1 3 

Eq. (4) can not be solved analytically because of the last term in the equation, we prefer 
to use a mathematical identity such as dm{r)/dr — —dT:{r)/dr to eliminate this term. We 
obtain the mass function from the identity as 

m{x) = mo + mie-^^ + mae-^'^^ , (31) 

where mo corresponds to the integral constant, and the parameters mi, and m2 are 2D, and 
—D, respectively. The parameter mo will denote the rest mass of the Dirac particle. By 
using the potential form given by Eq. (29) replaced by Eq. (28), inserting the mass function 
in Eq. (31), setting the "difference" potential A(r) to generalized Morse potential in Eq. 
(27) and using the new variable s = e~^^ , we have 

Id 1 

{t^ ^ ~T +~[~ '^^(«o«i + ml- E^) - 6^[aoa2 + (mo - E){mi + 2D)]s 
as s as s 

- 8^[aoa3 + (mo - E){m2 - D]s^] }(t>nn{s) = . (32) 

Comparing Eq. (32) with Eq. (13) gives the parameter set 

cci = 1 , -^1 = -(5^[aoa3 + (mo - E){m2 - D\ 

Q;2 = , ^2 = -5^[aoa2 + (mo - E){mi + 2D)] 

as = , -^3 = -5^(aoai + mg - E"^) 

-4 = 0, ^5 = ^^^^ 

OtQ^il, OL-J^ -^2 

as = 6 , ag = 

aio = 1 + 2VI3 , an = 2V|i 

ai2 = \/6, ai3 = -Vli 

where b — Ija.. We write the energy eigenvalue equation of the generalized Morse potential 
by using Eq. (26) 

ox/ 1 — i m . ^0(^2 + {mo- E){mi + 2D) 

2dyaoai +mQ — E^ — = = 2n + l . (34) 

'aoos + (mo - -E)(m2 - D) 
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Since the negative energy eigenstates exist in the case of the pseudospin symmetry [14, 15, 
16], so we choose the negative energy solutions in Eq. (46). In Table I, we give some nu- 
merical values of the negative bound state energies obtained from Eq. (46) for CO molecule 
in atomic units, where we use the input parameter set as = 11.2256 eV, Tq = 1.1283 A, 
mo = 6.8606719 amu, and a = 2.59441 [46], and summarize our results for different i, and 
n values. The corresponding lower spinor component can be written by using Eq. (25) 

0(s) = s-ie-^2.^fi(2w2s), (35) 

where Wi = S^J a^ai + ml — , and W2 = S^Ja^a^ + (mg — E){m2 — D) . 

Let us study the two special limits, pseudospin and spin symmetry cases, respectively, in 
the case of the constant mass. 



1. Pseudospin Case 

The Dirac equation has the exact pseudospin symmetry if the "sum" potential could 
satisfy the condition that dT,{r)/dr — 0, i.e. E(r) = A(const.) [14]. The parameters in our 
formalism become mi = m2 = 0. Setting the "difference" potential A(r) to the generalized 
Morse potential in Eq. (27), using Eq. (29) for the term proportional to 1/r^, and using 
the new variable s — e~^^, we have from Eq. (4) 

{ ^ ^ s ^ ^ ^ [ ~ ^ ^ ~ 5^{2MD + aoa2)s 

+ 5^{MD - aoaa)^^] }(/»(s) = . (36) 

where M = mo + A — E. By following the same procedure, the energy eigenvalue equation 
for the exact pseudospin symmetry in the case of constant mass is written 

/ : ~ and) + 2DM , , , ^ 

2 Jaoai + M mo + E) = J - + a{2n + 1 . 37 

and the corresponding wave functions read as 

^m^=m2=0 ^ e- ^'^^L'^ {2w'2s) , (38) 

where w'-^ — S^aoUi + M{mo + E) , and w'2 — 5\/aQa^ — DM . We must consideration the 
negative bound states solutions in Eq. (37) because there exist only the negative eigenvalues 
in the exact pseudospin symmetry [14, 15, 16]. 



2. Spin Case 



The spin symmetry appears in the Dirac equation if the condition is satisfied that A(r) = 
Vv{r) — Vs{r) — A{const.). In this case, we have from Eq. (3) 

- ^^^7^ -{mo + E-A){mo-E- E(r))}x(r) = , (39) 

where we set the "sum" potential as gcnerahzed Morse potential given in Eq. (27), and use 
approximation for the term proportional to in Eq. (29) [45] 

Vm{x) = boih + 626-^" + 636-'''") , (40) 

where 60 = '«(/« + l)/'^o > the parameters bi{i — 1, 2, 3) are given in Eq. (30). Using the 
variable s — e~^^, and inserting Eq. (40) into Eq. (39), we obtain 

- 5^(bob3 + DM')s^]}x(s)^0. (41) 

where M ' — mo + E — A. We write the energy eigenvalue equation, and corresponding wave 
equations in the spin symmetry limit, respectively, 
5[2DM'-bob2\ 



^/^_^_^ + 2VM, + M'(mo-^) =2n+l, (42) 
and 

X"''='^'=%s) = s< e- <-Lf 1 (2<s) , (43) 

where w'( = S^b^bi + M'{mo — E) , and w'2 = Sy/bob^ + DM' . We must take into account 
the positive energy solutions in Eq. (42) in the case of the exact spin symmetry [14, 15, 16]. 

In Summary, we have approximately solved the effective mass Dirac equation for the 
generalized Morse potential for arbitrary spin-orbit quantum number k in the position- 
dependent mass background. We have found the eigenvalue equation, and corresponding 
two-component spinors in terms of Lcgcndrc polynomials by using the parametric NU- 
method within the framework of an approximation to the term proportional to . We 
have also obtained the energy eigenvalue equations, and corresponding wave functions for 
exact pseudospin, and spin symmetry limits in the case of constant mass. We have observed 
that our analytical results in the case of the pseudospin symmetry are good agreement with 
the ones obtained in the literature. 
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TABLE I: Energy eigenvalues for the CO molecule for different values of i and (n, k) in the case 



of position dependent mass. 


e n 


K 


state 


< 


1 1 


-1 


lSl/2 


6.15913020 


2 1 


-2 


lP3/2 


6.52968379 


3 1 


-3 


lc?5/2 


6.89146288 


4 1 


-4 


1/7/2 


7.24974882 
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